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The decay rj —>■ 7r“''7r“7r° is studied with the KLOE detector, at the 
DAc^NE e“''e“ collider. Using a data sample corresponding to an integrated 
luminosity of 1.6 fb“^ a new study of the Dalitz plot is presented. 
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1. Introduction 

The isospin violating decay r/ —>• Svr is a good laboratory to study the 
light quark mass difference, — niu- According to Sutherland and Bell [U 
[2] , the electromagnetic contributions to this decay are small. Calculations of 
the electromagnetic corrections in chiral perturbation theory (yPT) at next- 
to-leading order (NLO), at the isospin limit, O(e^) [3] and including isospin 
breaking to leading order, 0{e^{md — mu)) [1], show only small effects, 
confirming Sutherland and Bell’s conlcusions. The 7 / —)• Svr decay thus 
proceeds mainly by strong interactions and the process can be calculated 
in yPT. Up to NLO the decay rate of r/ —)• vr+Tr^vr*^ is related to the quark 
mass ratio Q as follows: 


4_ _ n, ^-A ^9 ml — m^ ^ 1 , 

r(?7->7r’^7r 7r^)(xQ ^^ m =-(m^-b m„). 

m^ — m^ 2, 

As can be seen from the definition of Q, if one neglects the small term 
rh?/ml., the formula becomes an ellipse in the ^ plane with major semi¬ 
axis Q, the Leutwyler ellipse [5], shown in figure A precise knowledge of 
Q is an important constraint for the light quark masses. 

At leading order in yPT, Q can be extracted from meson masses. As¬ 
suming Dashen’s theorem |6] to fix the electromagnetic contribution to the 
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Fig. 1. The Leutwyler ellipse shows the contraints for the light quark masses, here 
for Q = 22.7 ±0.8. 


meson masses gives Q = Qd = 24.2. Using Qu in LO and NLO xPT 
calculations results in values of the decay width r(r/ —)• 7r'’'7r“7r*^) not in 
agreement with the experimental value: F^o = 66 eV, F nlq = 160 ± 50 eV 
[7] (updated to Ftvlo = 168 ± 50 eV [5]) and Texp = 295 ± 16 eV [8]. 

The big difference between LO and NLO results indicates a slow con¬ 
vergence of the xPT series and the importance of final state interaction 
between pions. Instead of using Q to predict the decay width, an accurate 
theoretical calculation of the decay amplitude together with an accurate 
experimental value for the decay width can be used to extract Q. 

There have been several efforts on improving the decay amplitude calcu¬ 
lations. A NNLO xPT calculation was performed by Bijnens and Ghorbani 
(in leading order isospin breaking) [^, but at this order many Low Energy 
Constants must be fixed. Dispersion relations can be used to take into ac¬ 
count the TTvr interactions, and two recent different methods have been used 
to improve the NLO xPT result mm- These promising methods can use 
as input the experimental Dalitz plot density information for the charged 
T] —)> 7r'*“7r“7r^ channel. Their predictions provide a precise determination 
of Q, and also provide an important consistency check, namely to get a 
reasonable agreement with the experimental results for the neutral channel 
r] —7- Svr^. As one of the main sources of uncertainty on Q comes from the 
uncertainty on the input parameters, it is crucial to precisely measure the 
density distribution of the r] —?• 7r'’'7r“7r^ Dalitz plot. 
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Experiment 

—a 

b 

d 

/ 

Gormley(70) [14j 
Layter(73) [T5] 
CBarrel(98) [16] 
KLOE(08) [121 
WASA(14) [I3| 

1.17(2) 

1.080(14) 

1.22(7) 

1.090(5)(;i9) 

1.144(18) 

0.21(3) 

0.03(3) 

0.22(11) 

0.124(6)(10) 

0.219(19)(47) 

0.06(4) 

0.05(3) 

0.06(fixed) 

0.057(6)(1[6) 

0.086(18)(15) 

0.14(1)(2) 

0.115(37) 


Table 1. Experimental results for Dalitz plot parameters of ry —>■ tt+tt 7r°. 


1.1. Dalitz plot for the charged decay 

Only few experiments since the 70s have measured the r/ —>■ 7r'’'7r“7r^ 
Dalitz plot distribution. The results for the extracted Dalitz plot parameters 
are shown in Table This Dalitz plot is usually presentend using the X 
and Y variables: 


X = v/3 
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Qrj — T ^7r~ T 2?7ryj-+ 


where T is the kinetic energy of the different pions in the p rest frame 
and s, t, u are the usual Mandelstam variables. The Dalitz plot parameters 
a,b,c,..., presented in the table, are defined by a polynomial expansion of 
the amplitude squared in X and Y, \A{X, Y)\^ ~ X(1 + oY + + cX + 

dX^ + eXY + fY^ + gX'^Y + hXY"^ + IX^) [) 

As can be seen in Table there is some tension between the two most 
recent, high statistics experiments (KLOE [12] and WASA m), and more 
data is needed to provide a clearer picture. A new KLOE analysis, using an 
independent data set, a better understanding of the detector and improved 
analysis methods aims to provide a new measurement with reduced system¬ 
atic uncertainties. The status and methodology of this analysis is presented 
here. 


2. Signal event selection 

The analysis is performed using 1.6 fb“^ of data collected with the KLOE 
detector in 2004-2005 (the previous analysis was based on 450 pb“^) . The 
KLOE detector is located in the DAijiNE (/)-factory, and the p meson is 


^ The parameters c,e,h and I would signal a violation of charge conjugation invariance 
and are therefore usually fixed to zero in the fit of the Dalitz plot distribution. 
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produced via the radiative decay (j) —)■ r]'y^. For the r] —>• Dalitz 

plot analysis, with —>■ 77 , we thus have two charged tracks and three 
photons in the final state. The event selection requires three clusters in 
the calorimeter m consistent with photon^ and two tracks with opposite 
curvature in the drift chamber |18j . 

To improve the signal to background ratio several cuts are applied: on 
the angle between the tracks and the photons; on the time-of-flight used for 
particle identification; on the angle between the decay photons in the tt^ 
rest frame (0jro); and on the missing mass MM{(j} — tt~^ — tt~ — 7 ,^). The 
two last cuts are depicted in Figure where both data and Monte Carlo 
simulation are shown. These histograms are used to fix the scaling factors 
used for the Monte Carlo background, and, as can be seen in Figurej^ a good 
data-Monte Carlo agreement is obtained. After reconstruction and cleanig 
cuts the signal efficiency is 37.6% with a 1% background contamination, 
evaluated from Monte Carlo. 


3. Fit of the Dalitz plot 

The Dalitz plot parameters are obtained by a fit of the parametrization 
shown in section o to the measured Dalitz plot density, by minimizing: 

Nbins / a/- _ ^Nbins ^ q j\t \ ^ 

^j=l ^j^ij^^theoryj ^ 


= 


2=1 


(Ti 


where Ni is the background subtracted data content in Dalitz plot bin i, ej 
is the acceptance of bin j, Sij is the smearing matrix from bin j to bin i, 
Ntheoryj = J \A{X,Y)\‘^dXdY (inside the phase space) is calculated with 
Monte Carlo integration for each bin, and ai, the error in bin i, includes the 
error in Ni and in ejSij. Figure]^ shows the experimental Dalitz plot used. 

From a detailed study both on data and Monte Carlo, fixing to zero the 
c and e parameters (and keeping the amplitude squared expansion up to 
the / term), it is found that the statistical uncertainty on the a,b,d and 
/ parameters can be reduced by about a factor of two with respect to the 
previous KLOE result (see Table [^, while improving also the systematic 
uncertainties. 

The inclusion in the fit of the g parameter would also be possible and 
its impact on the stability of the result is currently under study. 

In the near future, a further significant improvement of both statistical 
and systematic uncertainties of the Dalitz plot parameters will be possible 
at KLOE-2 [T9j, thanks to the luminosity upgrade of DA(^NE and to the 
better quality of reconstructed data using new detectors: inner tracker [ 20 ] , 

^ Not connected to a track and within the expected time window for a massless particle. 
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Fig. 2. Data-Monte Carlo comparison for the missing mass squared (top) and 
opening angle between tt^ photons in the 7r° rest frame (bottom). The vertical 
lines show the selection cuts. 


crystal calorimeters (CCALT) |21j and a tile calorimeter (QCALT) |22j of 
the upgraded KLOE detector. 
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Fig. 3. Experimental Dalitz plot of 77 —?► tt+tt tt®. 
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